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By Harold C. Lester 
Langley Research Center 


SUMMARY 


The stability characteristics of a representative launch vehicle utilizing 
a type of nonlinear thrust vectoring are investigated by the describing-function 
method. The nonlinearities in the system are associated with the control thrust 
forces which are produced by an auxiliary engine rotating through large gimbal 
angles with predetermined limits imposed on the amplitude of rotation. Limit 
amplitudes as large as ± 90 ° are considered and a describing function is derived 
to approximate the nonlinearities introduced by the amplitude limits and the 
sine coupling of the control engine with the rigid-body launch-vehicle equations 
of motion. Stability bounds for the system are determined graphically by means 
of a gain-phase analysis for various combinations of the control system and 
excitation parameters. Selected stable limit cycles, predicted by the 
describing function analysis, are verified by using an analog computer simula- 
tion of the launch- vehicle configuration. 


INTRODUCTION 


The control and guidance of almost all present-day launch vehicles is 
achieved by thrust vectoring, that is, by swiveling the thrust chambers of all 
or some of the main booster engines to produce the desired stabilizing force 
components. Maximum thrust-vector-rotation angles required are usually less 
than ±10°. Although such a control method has been successful for past and cur- 
rent applications, new thrust-vector control techniques may be more suitable 
with future launch vehicles, particularly the large solid-propellant systems. 

A recently proposed method employs auxiliary control engines producing a much 
lower thrust than the main booster engines. In order to obtain lateral forces 
and pitching moments large enough to stabilize and control the vehicle, the con- 
trol engines are rotated through much larger gimbal angles - as high as ± 90 °. 
Such a thrust vectoring scheme leads to nonlinear coupling of the control thrust 
with the vehicle equations of motion by virtue of the requirement for larger 
rotation angles. An additional nonlinearity is Introduced by amplitude limits 
which are imposed on the control engine. 

Linear stability methods - such as the Bode, Nyquist, and root-locus pro- 
cedures - are not directly applicable for the study of nonlinear systems. 



However, a linearization technique known as the describing-function method 
(ref. l) can be used in certain cases to approximate the response characteris- 
tics of a nonlinear system in order to take advantage of the well-established 
methods of linear analysis. This method has been applied with success in ana- 
lyzing servomechanisms involving such nonlinearities as off -on relays, satura- 
tion, deadband, coulomb friction, backlash, and hysteresis (refs. 1 , 2 , and 3 ). 
The method has also been used in analyzing off -on reaction control systems 
employed on the upper stages of vehicles for control and stability during the 
terminal phase of exit from the earth’s sensible atmosphere and for attitude 
control in orbit (ref. 4). In addition, the describing function has been 
employed to study the nonlinearities present in a particular element within an 
open- or closed-loop system (ref. 5 )» 

The present study examines the applicability of the describing-function 
method for investigating the stability of a launch vehicle having the aforesaid 
nonlinear thrust -vectoring characteristics. In the study, the gain-phase char- 
acteristics of the system nonlinearities are approximated by a describing func- 
tion. The remaining portion of the control loop which is linear is represented 
by gain-phase plots, and the stability of a typical launch-vehicle configuration 
is determined graphically. 

Motion in the pitch plane is considered and constant-coefficient differ- 
ential equations, which describe the perturbations about a discrete point in an 
ascent trajectory, are used to represent the rigid-body dynamics of the launch 
vehicle. The control engine is positioned in accordance with command signals 
from an attitude and attitude -rate feedback sensor in the control loop. Param- 
eters associated with a representative launch-vehicle configuration are evalu- 
ated at the maximum dynamic pressure condition and stability bounds are pre- 
sented for various combinations of attitude and attitude -rate gains, control 
engine frequency and damping parameters, and excitation parameters. Selected 
stable limit cycles, predicted analytically, are verified by using an analog 
computer simulation of the system. 


SYMBOLS 


A nondimensional amplitude of unlimited sinusoidal response of control 

engine ( eqs . ( 12 ) ) 

aQ control-system attitude gain 

a^ control-system attitude-rate gain, sec 

B nondimensional amplitude associated with sinusoidal response of con- 

trol engine with amplitude limits imposed (eq. (l4)) 

D drag force, lb 

g gravitational acceleration constant, in./sec^ 

2 



G(s) transfer function for linear elements in open-loop system (eq. (24)) 

I mass moment of inertia of launch vehicle about center of gravity, 

lb-sec 2 -in. 

k index of summation 

l length of launch vehicle, in. 

m total mass of launch vehicle, lb-sec 2 /in. 

N describing function (eq. (9)) 

S(x) local cross-sectional area of revolution of launch vehicle about 

X-axis, in. 2 

s Laplace transform variable, l/sec 

booster thrust, lb 

T c control thrust, lb 

t time, sec 

V velocity of launch vehicle, in. /sec 

X,Y body axis coordinate system 

x,y coordinates along X and Y axes, in. 

x coordinate locating center of gravity, in. 

y normal perturbation variable, in. 

(3 control thrust (control engine) rotation angle, radian 

3 C control thrust command function, radian 

P^ limit amplitude, radian 

3 q amplitude of sinusoidal forcing function (eq. ( 7 )), radian 

5 nondimens ional control engine rotation angle, 3/3 0 

7 = sin 3 

generalized coefficient associated with Fourier expansion of 7 
(eq. ( 8 )) 

0 attitude perturbation variable, radian 
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0 e attitude error angle, radian 

0J, feedback signal, radian 

0 phase angle associated with G(s), radian 

e> 

0 programed attitude, radian 

ir 

[i attitude-rate gain parameter, a o/ a l^ l/ sec 

£ viscous damping ratio associated with control engine 

p atmospheric density at altitude of interest, lb-sec^/in^ 

t nondimens ional time (eqs. (lO)) 

T n particular value of nondimensional time t for n = 0, 1, 2, 3> 

and 4 (see fig. (l 4 )) 

t 1 nondimensional time as measured from T2 j that is, t ! = t - T2 

a nondimensional amplitude ratio, & 0 /&i 

0 attitude associated with ideal ascent trajectory, radian 

0 n particular phase angle associated with sinusoidal response of control 

engine for n = 0, 1, and 2, radian 

i|r, generalized phase angle associated with Fourier expansion of 7 

* (eq. (8)) 

H nondimensional frequency ratio, ay cOq 

nondimensional damped frequency ratio, 03/03^ 

03 forcing function frequency (eq. (7))^ radian/sec 

03 n undamped natural frequency of control engine, radian/sec 

03 ^ damped natural frequency of control engine, radian/sec 

A dot over a variable indicates a differentiation with respect to time t. 
A primed variable indicates a differentiation with respect to x. Bars over 
symbols denote vectors. 
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ANALYSIS 


In this section the "basic equations of motion are presented and discussed. 
A describing function is then calculated to account for the system nonlinear- 
ities. Finally, the transfer function of the closed-loop system is derived and 
the governing stability equations established. 


System Equations 


Coordinate system . - Generally, 
mined by analyzing the perturbation 
discrete values of flight time. 

The situation is illustrated in 
figure 1 which shows the degrees 
of freedom employed in the present 
analysis: a normal translation of 

the center of gravity y; a rigid- 
body rotation about the center of 
gravity 0; and a control thrust 
rotation angle 0. It should be 
noted that the y and 0 degrees 
of freedom represent small devia- 
tions from an equilibrium flight 
condition. The resulting pertur- 
bation equations have time-fixed 
coefficients and are linear in the 
variables y and 0. 

Equations of motion . - The 
equations of motion for the trans- 
lation and pitch perturbations can 
be derived by the method of Lagrange 
(see, for example, ref. 6) and are 
summarized below: 


the stability of a launch vehicle is deter 
motion about a reference trajectory for 



Normal perturbation equation: 


y = 





(x) dx 


y + 


Tb - D 


m 


(x - x) S ? (x) dx 




S* (x) dx 


0 



P 


(i) 
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Pitch perturbation equation: 


0 = 


/ (x - x) S 1 (x) dx 

1 


y - 


f (x - x) S*(x) dx 
1 J 0 


“P r (x - x) 2 S'(x) dx 
1 J 0 


/? c x> \ 

0 + l— J-Isin 3 


( 2 ) 


In equations (l) and (2) the aerodynamic forces and moments have been 
developed by using an approximation of slender-body momentum theory (ref. 7). 
This theory gives the distributed aerodynamic lift (per unit length) as propor- 
tional to the derivative S T (x), where S(x) represents the local cross- 
sectional area of revolution of the launch vehicle about its longitudinal axis. 
One additional comment concerning the thrust forces is necessary. The total 
thrust component Ty acting in the y direction is given by the equation 

T- = (Th + T c cos p)e - T c sin p 

where T-^ and T c are the booster and control thrust, respectively, and 0, 
it should be recalled, represents a small perturbation in pitch attitude. With 
the restriction that the ratio d^/T^ be small relative to unity - for example, 
0.1 or less - the following approximation 

T- ~ T b e - T c sin p 

«7 

is assumed in deriving equation (l). 

Control engine . - Launch vehicles are controlled by lateral-force components 
which for the configuration considered are obtained from an auxiliary control 
engine. This method differs somewhat from the conventional thrust vectoring 
method. The two methods may be compared with the aid of figure 2. As illus- 
trated in part (a) of this figure, the conventional thrust vectoring method 
utilizes the thrust T-^ of the main booster engines by rotating the thrust 
chamber through a small angular misalinement or gimbal angle p. Maximum gimbal 
angles are usually less than ±10°. With the auxiliary control engine method, as 
shown in part (b) of figure 2, the main booster engines are fixed. Lateral 
thrust components are obtained in this case by rotating the control engine 
through gimbal angles which may be as large as ±90°. With the conventional 
method, the lateral thrust components are proportional to the gimbal angle 3 
(in radians) since rotation angles remain small. However, with the auxiliary 
engine method, the lateral thrust components are proportional to sin p and 
result in a nonlinear coupling with the rigid-body equations of motion. In 
addition, if the control engine is amplitude limited at some arbitrary level 
P i, where p ^ < ± 90 °, an additional nonlinear influence is possible. 
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1/3 l< io° 


|/3 1 ^ 90 ° 


(a) Conventional method. (b) Auxiliary control engine method. 

Figure 2.- Methods of thrust vector control. 


In the present analysis the control engine dynamics can he summarized hy 
the following second-order differential equation: 

3 + 210^3 + = ^p c (3) 

subject to the boundary conditions or constraint equations 

IMS' 3 ! ] 

9 = 0 (fi = ±9;) | (M 

where p^ is an arbitrary limit amplitude and p c is the command input to the 
control engine. (See appendix A.) It should be noted that for p < p£, the 
response of the control engine is that of a damped linear second-order system. 
For p = p^, however, the response is nonlinear as a result of the boundary 
conditions . 

Control-system considerations . - The control-system block diagram is shown 
in figure 3. As indicated in the figure, an attitude and attitude -rate feedback 
sensor is assumed so that in the time domain the total feedback 9f is given 
by the equation 
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( 5 ) 


0. = 0 + - 9 
f M- 


where \i is an attitude -rate gain parameter which is defined as the ratio of 
the attitude gain a q to the attitude -rate gain a^_, that is. 


M* = 


fl) 

a l 


(6) 


In the expression for the total feedback signal (eq. (5)) the 0 component 
serves to introduce damping into the rigid-body control mode. 


Nonlinear element, N 



Figure 3.- System block diagram. 

In addition, figure 3 illustrates the location of the system nonlinearities 
in relation to the closed-loop control-system structure. The nonlinearities are 
associated with the two forward path elements which are enclosed by the dashed 
line. The input to the total equivalent nonlinear element N is (3 C , the con- 
trol engine command signal. As indicated by equations ( 3 ) and (4), the output 
response (3 of the control engine may be either linear or nonlinear depending 
on whether the control engine limits ^(3 = ±|3 ^). Finally, the sin 3 coupling 
of the control thrust with the translation and pitch degrees of freedom is indi- 
cated in figure 3 "by the output signal 7 of the nonlinear element. In the 
ensuing section a describing function N is derived to approximate the 
frequency-response characteristics of the system nonlinearities. 


Describing-Function Procedure 

General procedure .- The Nyquist, Bode, and root-locus procedures which form 
the basis of feedback system analysis are limited in applicability to linear 
systems. Methods available for the analysis of nonlinear systems are not nearly 
as plentiful and as general in application as those developed for linear 
systems. The phase-plane method, for example, requires that a system be repre- 
sentable by a single second-order differential equation and only predicts the 
transient response from a given initial state. 

The difficulties encountered in nonlinear analysis are related to the 
exceedingly complex behavior of nonlinear systems as compared with linear 
systems. For example, the response of a linear system is linearly dependent 
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on the amplitude of applied excitation; doubling the input amplitude merely 
increases the amplitude of the output by the same multiple. With nonlinear 
systems, the output can be significantly dependent on the level or amplitude of 
excitation. Furthermore, if a linear system is excited by a sinusoidal forcing 
function, the output will be a sinusoid of different amplitude and phase but 
having the same frequency. The output of nonlinear systems forced in a similar 
manner may contain harmonics of the excitation frequency. One peculiarity com- 
mon to nonlinear systems which will be of interest in the present analysis will 
be constant -amplitude oscillations known as limit cycles. 

The describing-function technique is a linearization method based on the 
assumption that the input to the nonlinear element in a system is a sinusoid 
(refs. 1, 2, and 3 ). Graphically, this method is illustrated in figure 3 which 
shows the block diagram of the system under consideration. The nonlinear ele- 
ment shown here represents the totality of all nonlinearities contained within 
the system. The input to the nonlinear element is assumed as 

0 C = 0 Q sin cut ( 7 ) 

that is, a sinusoid of amplitude 0 O and frequency a>. The output y(t) is 
then expressed in the form of a Fourier sine series 


7 (t) = 7 1 s invest + ^ 7 k sin^teut + \|r^ ( 8 ) 

v ^ j 5' • • • j 

Fundamental Higher order terms 

where the coefficients 7 ^ (for k = 1,3*5* •••) are functions of both (3 Q 
and a). It should be noted that the even coefficients (k = 2,^,6,...) are 
identically zero since these terms evolve from an integration over the half 
period (0 < cut < rt) of the product of an even and odd function. The key assump- 
tion on which the validity of the describing function depends is that only the 
fundamental component of the output is fed back and contributes significantly 
to the input of the nonlinear element. (See ref. 2.) This assumption is per- 
missible on the grounds that the passive elements in the loop act as low-pass 
filters and attenuate the higher harmonics to a greater extent than the 
fundamental. 

The describing function N is defined as the ratio of the coefficient of 
the fundamental to the amplitude of the input (refs. 1 , 2 , and 3 ) : 


N(Po^) 


7l(Po> m ) /*l(P QJ 60 ) 

= 0 9 

F o 


(9) 


and can be interpreted as a linear element in the control loop which has an 
amplitude 0 Q as well as frequency-dependent cd gain and phase. When 7 ^ is 
independent of frequency, the describing function is said to represent a simple 
nonlinearity. In such cases, the phase angle of N is zero and only the gain 
of the system is affected. If the describing function is both amplitude and 
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frequency dependent., it is said to define a complex nonlinearity. The 
describing-function method enables the analyst to retain use of the linear tech- 
niques of analysis by extending these methods to compensate for the variable 
gain and phase of the nonlinear element. 

Derivation of the describing function.- Describing-function analysis of 
the control engine requires that the forcing input to the control engine be 
assumed as a sinusoid. Thus, it is necessary to determine a solution to equa- 
tions (5) and (4) with p c given by equation (7). By using the following non- 
dimensional parameters 


-(£> 




n = 


CJD 




n 


a = = 

a ov, 


n 


Vi - i 2 


O = 




(10) 


the sinusoidal response of the control engine is derived in appendix A and the 
solutions for both the linear and nonlinear response modes may be summarized as 
follows : 


Linear-response equations 



where 


S(t) = A sin^2rtr - 




A = 


^4| 2 fi 2 + (l - ft 2 ) 2 

0 2 = tan" 1 /'-^ 




1 - SI* 


Nonlinear-response equations 



(ID 


( 12 ) 


10 




I 



■-G) 


t a = ■t— sin 
0 2rt 


T, = T, + — 

5 12 


(a > 1 . 0 ) 


T 4 = t 2 + - 


1 

2 


(I 4 e) 

(I4f) 

(l^g) 


and where can he found hy solving the equation 

A 3 in(2«T 3 - 0 2 ) + ELF- ^i(T, - T 2 )ll .H.fe(Tj - t 2 ) + 0^1 + § - 0 (15a) 


subject to the conditions: 


|- T 0 

(cr > 1 . 0 ) (15b) 


(a < 1 . 0 ) (15c) 

given by equations ( 10 ) to (15), an 


analog computer solution of equations (5), ( 4 ), and (7) was effected. The 
results, for several input conditions, are presented in figure 4 . Part (a) of 
figure 4 illustrates the response obtained with a damping ratio of £ = 0.4 


when the system is excited near resonance 


0 ) 


=1.2 and forced with an ampli- 


tude which is twice the limit amplitude ~ = 2 . In part (b) of figure 4 , the 

° Z 

damping and frequency ratios have been decreased to £ = 0.2 and = 0.6 and 

vn 

the system excited by a level of excitation equal to the limit amplitude 

~ = 1 . In part (c) of figure 4 the resonant response of the system is illus- 
P Z 

trated for a damping ratio of | = 0.2 and a forcing amplitude equal to one- 

(3 

half the limit value — = 0.5. Since, in this latter case, the system is 

PZ 

excited at its resonant frequency, the system experiences limiting and the 
resulting response is cusped; that is, the control engine does not remain on 
the limit for a finite time. 


The preceding paragraphs have discussed the steady-state response of the 
control engine when excited by a sinusoidal forcing function. The describing 
function for the total nonlinear element (see fig. 5) must also account for the 
nonlinear coupling of the control engine with the translation and pitch degrees 
of freedom. Thus, the output from the total nonlinear element is given by the 
equation 
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7 - sin ^ 


(16) 


where 0 and 5 (eqs. (10)) are related by the following equation: 

0 - (0jCf)8 ( 17 ) 

The equations for 5 , 0 , and 7 (eqs. ( 10 ) to ( 17 )) were programed on a dig- 
ital computer and were evaluated for various combinations of the frequency ratio 

Q, = — , amplitude ratio a = — and for two values of damping £ = 0.4 and 
^ P l 

5 = 0.7. By using a digital integration routine, the amplitude 7 and phase 
angle ^ of the fundamental were computed and, as stated previously in equa- 
tion ( 9 )> the describing function for the system was determined. 


Results are presented in figure 5 which shows the variation of the ampli- 
tude of the inverse describing function l/N with its phase angle l80° - tq* 
The inverse form is used for convenience, as will be apparent in a later sec- 
tion. The results of figure 5 sire generalized in terms of the nondimens ional 
parameters a and Q (see eqs. (10)) representing an amplitude and frequency 
ratio, respectively. As is apparent from the figure, if the natural frequency 
of the control engine is specified and a particular limit amplitude value 
(3^ chosen, the describing function or rather its inverse can be considered a 

function of the excitation amplitude and frequency, that is, a function of 0 Q 
and co. It should be noted that the describing function possesses both an 
amplitude ratio and a phase angle and therefore alters both the gain and phase 
of the system. In addition, as a result of the amplitude limits, the response 
of the fundamental possesses a jump region in the gain-phase plane. 



(a) p* = * ^ ^ = 0 * 4 . 

Figure 4.- Response to sinusoidal input of control engine with an amplitude limit. 
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w n 


Figured- Concluded. 





Derivation of System Transfer Function 

In a purely linear control system, the transfer functions of the individual 
elements may he cascaded in any order to produce the net gain and phase shift 
between two points in the control loop. Similarly, interpreting the describing 
function as an equivalent linear element with an amplitude- and frequency- 
dependent gain and phase leads readily to the formulation of the open-loop 
transfer function of the system analyzed herein. 

From figure 3 the transfer function for the forward loop may be written 
as follows: 


9(s) _ y(s) 9 (b) 

9 e (s) 9 e (s) r(s) 


(18) 


By substituting aQN for 7(s)/0 e (s) and letting N the describing function, 
represent the combined gain-phase characteristics of the nonlinear element, the 
preceding equation for 0(s)/0 e (s) becomes 


9(a) 

9 e (s) 


aQN 


9(b) 

7(s) 


(19) 


The transfer function 0(s)/7(s) can be obtained from the rigid-body 
equations (eqs. (l) and (2)) by substituting the results of equation (l 6). 

If the Laplace transform is taken with respect to the variable t, all initial 
conditions assumed to be identically zero, and the terms collected, equa- 
tions (l) and (2) can be arranged into the following matrix form: 


a ll a 12 
&21 a 22 


y(s)l j a 13 | 

> = 7(b) < 
0(s)J | a 23| 


( 20 ) 


where the coefficients a^j (where i and j are 1, 2, . . . ) may be deter- 
mined by carrying out the transformation and identifying like terms in equa- 
tions (l), (2), and (20). If the determinants A(s) and C(s) are intro- 
duced and defined as 


A(s ) 


C(s) 


a ll 

a l3 

a 21 

a 23 

a ll 

a 12 

a 21 

a 22 
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The required transfer function may be written as follows : 


e(s) = a(s) 
7 (s) C(s) 


( 21 ) 


The transfer function for the feedback path is given by the following 
equation: 

0 _p ( s ) i 

-7 -r- = 7f(s + n) (22) 

0(s) M. 

and hence the open-loop transfer function of the system is simply 


9f( s ) _ 0 ( s ) 9f( s ) 

e e (s) 0 e (s) 0(s ) 

which becomes on substitution of the results of equations ( 19 ) and. (22) 

9 f( s ) = A(s) / n 

e e (s) " n c(s)' 


By letting 


G(s) = -p ^|(s + n) 
H C(s) 


(23) 


(2k) 


and substituting s = jto, where j = the open-loop transfer function may 

be written in a more usable form as follows: 


e f (<n) 




(25) 


It should be noted that G(s) represents the transfer function of the linear 
elements in the open-loop system. (See fig. 3 - ) In a similar manner, N 
represents the equivalent transfer function for all the nonlinearities in the 
open loop. 


Method for Determining System Stability Graphically 

By using equations (19) and ( 25 ), the closed-loop transfer function 
0(s)/0 p (s) for the system may be formulated as 

9(b) _ 9(s)/9 e (s) 

0p( s ) 1 + N G(s) 


( 26 ) 


IT 


where the numerator is the transfer function of the elements in the forward 
loop. The characteristic stability equation for the system is obtained by 
equating the denominator of equation (26) to zero; that is, 

1 + N Gr(s) = 0 (27) 


Solutions of this equation represent conditions of marginal (neutral) stability. 
Rearranging equation (27) and expressing the results in a polar form produces 
the following equivalent equations which may be solved graphically: 


|g(co)I = 1 , 

|n( P o ^)| 

e g (a 3 ) - 180 0 - 


(28) 


where the substitution s = jca has been made. 


E 

< 



Phase angle, 


{a) Typical gain-phase plot for GM. 



1 1 1 1 1 

Phase angle, 180°— 1 // ( 


(b) Inverse describing function l/N. 

Figure 6.- Illustration of graphical solution in 
gain-phase plane. 


In order to effect a graphical solution 
of equations (28), the amplitude |G| is plotted 
against the phase angle 0g with the frequency 
a) as a parameter. A typical plot is illus- 
trated in part (a) of figure 6 for several 
values C^, and of the attitude -rate 

gain a^. The vector is drawn tangent to 

the = C 2 curve and is directed in the 
sense of increasing frequency. An inverse 
describing-function plot is then dimens ion- 
alized by selecting a value for the natural 
frequency of the control engine and by 

using equations (10) to convert the constant 
curves to constant 00 curves. A typical 
inverse describing-function plot is shown in 
part (b) of figure 6 for several values of the 
amplitude ratio cr^, cr^, .... The vector 

V 2 is drawn tangent to the 00 = curve and 
is directed in the sense of increasing ampli- 
tude ratio or ^assuming o\± > cr^ > 0*2 ^ °l) 
where a is defined by equations (10). 

Once the inverse describing function has 
been dimensional! zed, equations (28) may be 
solved graphically (ref. 3) by superimposing 
parts (a) and (b) of figure 6. Solutions are 
represented by those intersections which occur 
at a common frequency a). The procedure is 
illustrated in part (c) of figure 6. An inter- 
section, that is, a solution, is indicated for 
a frequency of 00 = 0 ^. The amplitude ratio a 
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for the solution is between cr = and 

a - aj and the attitude-rate gain is 

a l = C 2* 

c 

Solutions of equations (28) determine A 

conditions of marginal stability known as 5 : 

limit cycles. If a limit cycle is stable, J 

the system will return or converge, after a | 

slight disturbance, to the marginal stability 
condition. If the system diverges, the limit 
cycle is unstable. 


°r c 2 



Phose angle, 0g (u) ond 180° -'/ / j 


In general, a rigorous definition of sta- 
bility, that is, whether a solution of equa- 
tions (28) represents a stable or unstable 
limit cycle, is difficult to formulate. 
Instead, an argument is presented based on an 
intuitive extension of the definition of sta- 
bility for linear systems. The argument pro- 
ceeds as follows: Let G(od) represent the 

open-loop transfer function of a linear sys- 
tem. By using the Nyquist stability method, 
Gt(gd) is plotted in the complex plane (for 
0 < jcu < joo) as shown, for example, in 
part (a) of figure 7« A circle of unit 
radius is then constructed about the origin 
and the intersection of this circle with the 
negative real axis defines the critical point 
-1 + JO. Within the framework of the abbre- 
viated Nyquist method (if it is assumed that 
the open-loop transfer function G (o>) is 
stable), the relative stability of the closed- 
loop system may be determined by the amount 
of phase margin \|/ present when G(cju) = 1 
or rather, when G(u>) intersects the unit 
circle. For example, at point A in part (a) 
of figure 7 the phase margin is positive 

(by definition) and the system is stable, 
since the Nyquist plot of G(co) does not 
encircle the critical point -1 + JO. If 
G(cd) intersected the unit circle at point B 
where the phase margin \|r-g is negative, the 
closed-loop system would be unstable. In 
summary the system is stable if the phase 
margin is positive and unstable if the phase 
margin is negative. The condition of mar- 
ginal or neutral stability occurs for zero 
phase margin. 


(c) Typical graphical solution resulting when the 
gain-phase plot and the inverse describing 
function are superimposed. 

Figure 6.- Concluded. 


Imaginary 
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<b) Illustration of the describing function 
method using a Nyquist plot. 

Figure 7.- Comparison of the Nyquist stability 
method as applied to linear and nonlinear 
systems. 


The preceding ideas may be extended 
intuitively to a nonlinear system with an 
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amplitude- and frequency-dependent describing function. In this case, the crit 
ical point is located in the complex plane by the critical vector 

X /l80° - ♦ I (p OJ m) 

|N(e o ,o»)|' 

where it is apparent that both the magnitude and direction of this vector vary 
with and cd. As illustrated in part (b) of figure 7 , a solution is shown 

at point A. Since the solution represents a point of marginal stability (zero 
phase margin) the critical vector and the G(cd) vector coincide (G is the vec 
tor which traces out the G(co) locus). 

Suppose that position A is a solution and that the describing function 
N(0 o ,cd) is such that a minute increase in amplitude moves the critical 

vector to position B. Based on the ideas presented previously for linear 
systems, one is tempted to interpret the angle as a small increment of 

positive (stable) phase margin and to conclude that the critical vector has 
shifted its direction so that the closed-loop system is stable for increasing 
amplitudes. If the system is stable, the solution will eventually return or 
converge back to position A and, on this basis, the solution may be classified 
as a stable limit cycle. Conversely, with decreasing amplitudes, the critical 
vector will shift to position C where the phase margin is negative and the 

system is unstable. As a result, the system will diverge and with increasing 
amplitudes the solution at position A will again be reached. For a solution to 
be classified as an unstable limit cycle, it appears that increasing 
(decreasing) amplitudes must produce negative (positive) phase margins. 

Although the preceding argument lacks a rigorous proof, it appears to par- 
allel similar types of reasoning given in the literature (refs. 1, 2, and 3) 
and is probably applicable to a wide variety of situations. In any case, as an 
analog computer simulation confirmed, this type of argument did prove to be 
adequate in classifying the stability boundaries for the system and nonlinear- 
ities discussed in this paper. 

If reference is made to part (c) of figure 6 and this type of argument is 
used, the solution (intersection) may be classified as either a stable or 
unstable limit cycle. By assuming that Pj > a 2, it is apparent that a small 
increase in the_amplitude ratio p moves the critical vector in the direction 
of the vector Vp and results in a small increment of positive phase margin 
since 0g is now slightly greater than the phase angle l80° - ® ie solu- 

tion is thus a stable limit cycle. 

To conclude this section, the following rule (from ref. 3 ) > which may be 
helpful in classifying limit cycles, is given. The rule is paraphrased here to 
fit the specific situation illustrated in part (c) of figure 6. Consider an 
observer on a curve of constant attitude-rate gain - in this case the a^_ = C2 
curve - facing in the direction of increasing frequency, that is, in the direc- 
tion indicated by vector V3. ^The vector V^, it should be recalled, is 


20 



tangent to the = Cj_ curve. j Now consider the vector V 2 which originates 
on the inverse describing-function plot and is tangent to a line of constant 
frequency - in this case the cjd = curve. (it should he recalled that V 2 
is directed in the sense of increasing amplitude a.) The_rule may he stated 
as follows : If, to an observer looking along the vector V-^ the vector V 2 

crosses from left to right, the solution (intersection) represents a stable 
limit cycle; if it crosses from right to left, the limit cycle is unstable. 


DISCUSSION OF RESULTS 


In this section the describing-function analysis outlined in the preceding 
section is employed to determine the stability characteristic of a hypothetical 
launch vehicle for the maximum dynamic pressure condition of a nominal ascent 
trajectory. Stability bounds are established for small control engine rotation 
angles, when the response of the system is essentially linear, by the root-locus 
method. Then by using the inverse describing function computed previously and 
representing the linear portion of the open-loop transfer function by gain-phase 
plots, stability bounds for the nonlinear mode are determined graphically. 
Finally, an analog computer simulation of the launch -vehicle system is used to 
verify the describing -function results. 


Physical System 

The hypothetical launch vehicle analyzed in this study consists of a multi- 
stage configuration with an overall length to first -stage diameter ratio of 

about 8 . 3 . The thrust-weight ratio is about 1.25 and the control thrust level 

T 

was selected as 10 percent of the total booster thrust, that is, — = 0.1. The 

T b 

control engine is allowed to rotate through gimbal angles as large as ± 90 °, but 
may be arbitrarily limited at any amplitude up to and including this level. 
Values for the coefficients required in equations (l) and (2) are given in 
appendix B. 


Linear Stability 

For small amplitudes of oscillation the response of the control engine is 
linear and is governed solely by equation ( 3 )- In addition, the coupling of the 
control thrust with equations (l) and ( 2 ) may be satisfactorily approximated by 
replacing sin 0 with 0 (in radians). The transfer function for the control 
engine is therefore given by the equation 


B(s) = 

Pc( S ) s 2 + 2 f;CD rl S + Oi 2 
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Athtude-rate-gain parameter, ^ Attitude- rate-gom parameter, 


which may he cascaded with the remaining elements of the control loop (eq. (2k)) 
to produce the open-loop transfer function of the linear system 


6f(s) 

_ a cP§ 

A(s) (s + n) 

e e (s)_ 

Linear 

C(s) (s 2 + 2g<%S + cuj 2 ) 


The closed-loop linear stability of the launch-vehicle configuration was 
determined by the root-locus method (ref. 2). The locus of roots is determined 

by equating the denominator of the 
closed-loop transfer function to 
zero, that is. 


5 r 


(a) 


00 . 7 . 



(a) 5 = 0.7. 



(b) C = 0.4. 

Figure 8.- Linear stability boundaries. 


9f*( s ) 


0 e( s )-! Linear 


1 = 0 


and plotting the roots as a function 
of the open-loop gain in the complex 
plane . 

By using the root -locus method, 
it is possible to establish regions 
of stability in terms of the control 
parameters slq and \±. These 
results are presented in figure 8 
for several engine frequencies o^. 
Part (a) of figure 8 presents the 
stability boundaries for an engine 
damping of | = 0-7* With an atti- 
tude gain of a 0 = 10 and an engine 
frequency of = 2, the vehicle is 
stable for 0.35 1.2. If the' 

engine frequency is increased to 
= by the limits on \i become 
about 0.13 < M- ^ 2.7. For a given 
engine frequency, for example, 

= 3j th e attitude gain may only 
be increased to about a q = 3^ at 
which the vehicle is unstable for 
any choice of the attitude-rate gain 
parameter 1 1 . Hence, attitude gains 
greater than a q = 3 ^ require 
engine frequencies greater than 
ODn = 3-0. The results presented in 
part (b) of figure 8 for | = 0.4 
may be interpreted in a similar man- 
ner. The choice of \i for a 
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particular aQ determines the damping of the pitch mode, the smaller values 
of p producing better damped responses. It should be recalled that the atti- 
tude rate gain a^ is inversely proportional to p as indicated by 
equation (6). 


Nonlinear Stability 

In general, disturbances such as winds and maneuver requirements will 
necessitate larger gimbal angles. It is therefore necessary to ascertain the 
stability characteristics of the system for the nonlinear mode when larger rota- 
tion angles are required. Since the describing function, as illustrated in fig- 
ure 5, possesses both a gain and a phase angle, the root-locus method is too 
cumbersome for efficient use. Instead, a graphical solution is effected in the 
gain-phase plane as has been previously outlined in the " Analysis" section. 

The following numerical example is presented as being illustrative of the graph- 
ical method for determining stability bounds of the nonlinear mode. 


A typical variation of gain with 
phase for G(<n) is presented in 
figure 9 for an attitude gain of 
aQ = 20 . By using the results of fig- 
ures 5 and a solution of equa- 
tions (28) can be briefly demonstrated 
(see discussion following eqs. (28)) 
by considering, for example, the fol- 
lowing choice of parameters: 

a) = 2 
n 

£ = 0.4 

a Q = 20 

Since an engine frequency has been spec- 
ified, the inverse describing-function 
plot illustrated in part (b) of figure 5 
may be dimens ionalized in terms of the 
frequency cjq by using equations (10). 
Hence, constant 0 curves are converted 
to constant o> curves. A graphical 
solution can now be obtained by super- 
imposing part (b) of figure 5 on fig- 
ure 9 and determining the intersections 
that occur at a common frequency. 

One such Intersection is illus- 
trated in part (a) of figure 10 for an 
attitude -rate gain of about = 35. 

The corresponding frequency is cd = 2 



Figure 9.- Typical gain-phase plot of G(w) for aQ = 20. 
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Phase angle, deg 


(a) Stable limit cycle. 



(b) Unstable limit cycle. 

Figure 10.- Typical limit cycles illustrating the graphical 
solution for a Q = 20, w = 2, and £ = 0A 


and the amplitude ratio is about a = 0.9* 
The nature of the intersection may be 
determined by applying the rule given pre- 
viously. Recall first that the vector 
originates on the G(a)) plot and is con- 
structed tangent to the a^_ = 35 curve and 


points in the direction of increasing fre- 
quency ay, whereas the vector V 2 origi- 
nates on the l/N plot and is drawn tan- 
gent to the (jo = 2 curve and points in the 
direction of increasing amplitude a. 

Since, to an observer looking along the 
V 1 vector, the vector V 2 crosses from 
left to right, the intersection (solution) 
is a stable limit cycle. Therefore, when 
the system is operating with this combina- 
tion of parameters, that is, a Q = 20, 
a^_ = 35 j oiji = 2, and £ = 0.4, a steady- 
state oscillation of amplitude ratio 

P Q 

— = 0.9 and frequency oj = 2 radians per 

second will develop and will persist indef- 
initely. Note should be made of the fact 
that the amplitude p Q of the input signal 
P c to the nonlinear element (fig. 3) can 
only be determined once a finite value has 
been assigned to the limit amplitude 3^; 
for example, if the control engine is 
limited at (3^ = ±45°, then the limit cycle 
is such that the input amplitude is about 

(3 = ^0.5°. 


For the particular choice of input 
parameters, intersections at other fre- 
quencies are possible and, indeed, by 

finding the intersections occurring at many discrete frequencies, it is possible 
to generate a stability boundary. For example, in part (b) of figure 10 another 
solution is indicated for an attitude-rate gain of about a^ = 12 with a = 1.5 
and co - 0.8. Here the V 2 vector crosses the vector from right to left 

and indicates, in accordance with the rule, that the solution is an unstable 
limit cycle. Thus, for amplitudes such that < 1-5* the system is stable 

whereas for P> Q /$i >1-5* the system is unstable. 

Stability boundaries for the nonlinear response mode may, in general, be 
categorized in accordance with the four regions defined in figure 11. As noted, 
the regions are illustrated for a particular value of slq and are defined in 
terms of the attitude-rate gain aj_ and the amplitude ratio a which is 
defined as the ratio of the forcing amplitude p Q to the limit amplitude 0^ 
(eqs. (10)). Along the cr = 0 axis, the stable region is bounded on the left 
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by a lower limit on ap. For attitude -rate 
gains less than this value, the system is 
unstable. On the right, that portion of 
the boundary falling along the a = 0 axis 
is bounded by an upper limit on ap above 
which stable limit cycles occur. It should 
be noted that this portion of the stability 
boundary, that is, along a = 0, can be 
determined from a linear analysis and is 
comparable to the information presented in 
figure 8. In addition, an upper bound is 
imposed on the region of stability by the 
occurrence of unstable limit cycles. This 
upper bound manifests itself in the form of 
an amplitude ratio (a) level for which the system becomes unstable if exceeded. 
For amplitude ratios less than this critical level, the system is stable. 

If the region of stability does not intersect the a = 0 axis, as shown 
by the dashed line in figure 11, the system response is either unstable or char- 
acterized by a stable limit cycle and is indicative of a completely undesirable 
set of control-system-parameter values. Control parameters should be chosen so 
as to preclude the occurrence of stable limit cycles in order to insure that the 
system will converge back to the equilibrium, represented by the cr = 0 axis, 
when the disturbing influences terminate. 



Figure 11.- Classification of stability boundaries for 
a constant a Q . 


Results are presented in fig- 
ure 12 for other solutions of equa- 
tions ( 28 ) when other values of the 
parameters under consideration are 
scrutinized for possible limit cycles. 
Here a damping ratio of £ = O.t is 
considered and the amplitude ratio a 
is plotted against the rate gain ap 
with the frequency crossplotted 
as a parameter. Graphical solutions 
are indicated by the solid curves for 
four engine frequencies c% = 2, 3> 4, 
and 5 and for several different val- 
ues of the attitude gain aQ. 

Part (a) of figure 12, which 
shows solutions for £ =0.4 and 
aQ = 10, may be interpreted as fol- 
lows: Consider, for example, the 

curve for = 2. With ap = 8 the 
amplitude ratio is about a = 2, the 
highest value considered. Analysis 
in the gain-phase plane reveals this 
point to be an unstable limit cycle. 
Hence, the system is divergent for 
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(a) a Q = 10. 



(b) a Q = 20. 

Figure 12.- Nonlinear stability boundaries for = 0.4. 
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a > 2 and convergent for a < 2. The corresponding operating frequency is 
about cd = 0.5 radian per second. Decreasing the rate gain ap produces a 
decrease in the amplitude ratio until finally at about an ap = 4.7 and. 
a = 0.95^ the curve drops almost vertically to a = 0. This lower limiting 
value of the rate gain (on a = 0 ) represents the limiting point between the 
stable and unstable modes of the linear system (fig. 8 ) for aQ = 10 , = 2 , 

and g = 0.4. Further decrease in the rate gain ap requires an increase in 
the engine frequency cjc^ for stability of the linear system to be maintained, 
as was ascertained previously. Although not shown in part (a) of figure 12 , it 
should be noted that at a higher level of the attitude-rate gain ap, the 

cd^ = 2 curve will break away from the a = 0 axis as will the = 3 , 4, 

and 5 curves. (See, for example, parts (c), (d), and (e) of fig. 12.) The 
value of rate gain at which this breakaway occurs is equal to the upper limit 
value of ap determined from the linear analysis. 

In part (b) of figure 12, the effect of increasing the attitude gain aQ 

may be noted. Here the = 2 curve is again multivalued with the amplitude 

ratio or; and instead of approaching a discrete value of at cr = 0 , it 

reaches a lower limiting value of about ap = 11.4 and then reverses direction, 
a remaining almost constant with increasing values of a^. Whereas the upper 
part of the curve denotes an unstable (divergent) limit cycle, the lower part 
indicates a stable (converging) limit cycle. Increasing the engine frequency 
to a) = 3 produces a curve of similar shape, but which intersects the a = 0 

axis. The lower limiting value of rate gain 
is about ap - 6.3 which agrees well with 
the linear results (fig. 8 ). The curves for 
(% = 4 and = 5 also give lower limiting 
values for ap which agree well with the 
values obtained from the linear analysis. 

Stability boundaries for a 0 = 30 , 40, 
and 50 are shown in the remaining parts of 
figure 12. For example, part (c) of fig- 
ure 12 , which shows solutions for an attitude 
gain of aQ = 30 , indicates that a stable 
limit cycle would dominate the response for 
an engine frequency of = 2 and attitude- 
rate gains greater than 20. The linear 

analysis indicates the system to be unstable 
for this combination of parameters. Hence, 
the response is divergent for small ampli- 
tudes, but as the amplitudes grow, the 
response converges to a constant-amplitude 
oscillation with an amplitude ratio of about 
a = 1, depending on the rate gain. If the 
engine frequency is increased ^ = 3 j 

the lower limit on the rate gain ib about 
ap = 11 (along cr = 0). For rate gains 




(d) a Q = 40. 
Figure 12.- Continued. 
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lower than about a^ = 11, the system is 
unstable. The upper limit for stability is 
about a^ = 2 6 and larger rate gains produce 
stable limit cycles. It should be noted that 
in the latter part of the region, the con- 
stant CD curves (dashed lines) approached 
the constant c% curves (solid lines) almost 
tangentially. For this reason the frequency 
cd is not crossplotted in this area and the 
curves are faired. Increasing the engine 
frequency to ^ and 5, respectively, 

gives even greater latitude in the choice of 
the rate gain a^. Of course, an upper limit 
is imposed on the amplitude cr which, if 
exceeded, throws the time-fixed system into a 
diverging (unstable) mode. With attitude gains of a q = k-0 and 50 as illus- 
trated in parts (d) and (e), higher engine frequencies, for example cn n = 4, 
are required to preclude the occurrence of stable limit cycles. 

Increasing the damping ratio to | = 0.7, as summarized in figure 13, has 
the effect of shifting the curves to the right so that the stable limit cycles 
occur at higher values of the attitude-rate gain a^_, although it appears that 



(e) a Q = 50. 
Figure 12.- Concluded. 



Figure 13.- Nonlinear stability boundaries for £ = 0.7. 
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slightly higher engine frequencies are 

required for linear stability. Near the 
lower limiting values of the rate gain a^, 
where a transition from an unstable to a sta- 
ble limit cycle is possible, the gain-phase 
analysis is inaccurate as to whether an inter- 
section represents a converging or diverging 
mode. However, increasing the rate gain in 
order to move away from these marginal cases 
produces more discernible intersections. At 
these higher values of a^, an analog-computer 
simulation of the system provided a check on 
the stable limit cycles as given by the gain- 
phase analysis. Analog computer checkpoints 
are indicated by circles and are discussed 
briefly in the following section. 


Analog-Computer Verification 

An analog -computer simulation of the launch vehicle was used to check the 
stability bounds obtained by the describing-function analysis. Analog-computer 
checkpoints are shown in figure 12 by circles and indicate stable limit cycles. 
For example, in part (d) of figure 12, which shows the resulting stability 
boundaries for a Q = ^0 and g = 0.4, there is a checkpoint at about a^_ = 52 
and an amplitude of about a = 0.95* This point is associated with an engine 
frequency of = 5 and indicates a cyclic frequency for the limit cycle of 
about cjd = 5 radians per second. As is apparent from the figure, the agree- 
ment is very good. Decreasing the rate gain, there are several additional 
checkpoints at frequencies of a) = 2.5, 2.0, and 1.8 radians per second for an 
engine frequency of (% = 3 and the agreement is again good. Further decrease 
in the rate gain a-^ did not yield as good agreement since these points fall 
closer to the condition of marginal stability. In addition, it was not possi- 
ble to check the upper branch of these curves, which represents unstable limit 
cycles, because the effective damping is very low along this part of the 
boundary. 


(e) o 0 * 50 


S 

? if- 



Attitude rote gain, a. 


(e) a 0 = 50. 
Figure 13.- Concluded. 


CONCLUDING REMARKS 


A stability analysis has been conducted on a launch vehicle having a type 
of nonlinear thrust vectoring. The frequency-response characteristics of non- 
linearities, which arise from control forces produced by an amplitude-limited 
engine rotating through large gimbal angles, were approximated by a describing 
function. By representing the linear portion of the closed-loop system by 
gain-phase plots, stability bounds were determined graphically for a range of 
parameter values. Selected stable limit cycles, predicted by the describing 
function analysis, were verified by simulating some of these sustained, 
constant-amplitude oscillations on an analog computer. The frequency and 
amplitude of these stable limit cycles were in good agreement with the values 
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measured from the analog time histories. The small amplitude results were 
checked hy linearizing the equations and doing an independent root-locus anal- 
ysis. An interesting aspect of the study was that the describing function for 
the nonlinearities investigated exhibited a type of jump behavior. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 

Va. , June 22 , 1965* 
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APPENDIX A 


SINUSOIDAL RESPONSE OF CONTROL ENGINE 


This appendix is devoted to determining the sinusoidal response character- 
istics of an amplitude -limited-control engine. The response equations are 
derived for both the linear and nonlinear modes. 


Analysis 


For the analysis presented herein, the characteristics of the control 
engine dynamics are given by the second-order constant -coefficient linear dif- 
ferential equation 

0 + 2^(3 + = o^p o sin cut (Al) 

subject to an amplitude constraint (boundary condition) as defined by the fol- 
lowing two equations: 

|P| < 0 Z (A2) 

0 = 0 (0 = ±0 Z ) (A3) 


Equation (Al) can be interpreted as describing the dynamical state of a 
linear spring-mass oscillator with an attached viscous dashpot (damper) when 
excited by a sinusoidal forcing function with amplitude (3 0 and circular fre- 
quency a). The parameter may therefore be thought of as the undamped 

natural frequency of the system, £ the viscous damping ratio, and p^ the 
limit amplitude. For certain combinations of the input parameters, that is, 

|, and the unlimited linear response may exceed a particular 

value of under consideration. Solution of equation (Al) requires, under 

these circumstances, a joint consideration of the boundary conditions 
(eqs. (A2) and (A3)) and the resulting response of the system is nonlinear. 

If, however, |p| ^ P^, the system behaves linearly. 

In a physically realizable system a dissipative mechanism (buffer) would 
be employed with the amplitude limits. The effect of any such buffer arrange- 
ment has been neglected herein with the assumption that the velocity (3 is 
instantaneously reduced to zero when limiting occurs (when p =±p^). As will 
be apparent later, this assumption produces an output waveform which is piece - 
wise continuous for the nonlinear response mode. 

The simultaneous solution of equations (Al), (A2), and (A3) yields the 
desired response for any combinations of the damping ratio £, frequency ratio 
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00 / 0 ^, and amplitude ratio 3 o /0j* The solution may be generalized by substi- 
tuting the following nondimens ional parameters: 



Differentiating 3 = 3 Q 6 and using r = ^^jt produces 

d|3 — — = ££ 

dt ° dr dt 2jt dT 




(A5a) 


and a second differentiation produces, in a similar manner, the following: 


d 2 3 = B /jr_\ 2 d 2 6 


dt £ 


>V 2 * 


(A5b) 


dT c 


Equations (Al), (A2), and (A 3 ) may therefore be rewritten in a nondimens ional 
form by using the results of equations (A4) and (A5). 



+ 26 



+ 
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(a6) 

(AT) 

(A8) 
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APPENDIX A 


Unless stated otherwise, equations (A 6), (A7), and (a 8) will hereafter he 
referred to as "the system" where it is assumed that nondimensionality is 
implied. 

The solution of the system equations may he conveniently divided into two 
parts governing the linear and nonlinear modes of operations. The linear 

response is determined first. If 5 < ^ < (3^, the constraint equa- 

tions (A7) and (A8) do not influence the solution. Hence, equation (A 6) may he 
solved hy using the methods of linear ordinary differential equations or opera- 
tional techniques (Laplace transforms) and is given hy the following equations: 

5(r) = A sin(2rtr - 0 2 ) (b < j) (A9) 

where 

A . 1 (A10) 

lfin 2 fl 2 + (1 - o 2 ) 2 


0p = tan 


-i ( 26 n 

li - of 


(All) 


as a piecewise continuous function. 


Equation (A9) is the steady-state linear response of the system when excited 
with a sinusoidal forcing function and is applicable only so long as A ^ l/cr. 

It will now he shown that the nonlinear case (A > l/cr) can he constructed 

Figure 1^ illustrates a typical output 
waveform when the system is operating in 
the nonlinear mode. The derivation begins 
hy noting that as soon as the system 
leaves the limit at t ~ and r = t^_, 

it is once again linear. Thus, this por- 
tion of the solution can he obtained hy 
the methods of linear analysis. 



Consider a shifted time reference 
noted as r T (where t' = t - T 2) in fig- 
ure lk. At t t = 0, the state of the 
system hy virtue of the amplitude con- 
straint is given as follows : 




Figure 14.- Typical steady-state output waveform of the 
amplitude-limited control engine. 


5 = i 

cr 


0 


(A12) 


dT 


j 


32 



APPENDIX A 


As the system leaves the limit at t t = 0, the fall-off response can he deter- 
mined hy solving the differential equation 


d% /2>t\dS 

^72 + 2 H~jd7 


- + (tP “ (f t sln [ 2lt ( T ' + t 4 


subject to the initial conditions expressed hy equations (A12). The general 
solution may he formulated hy the straightforward methods of linear analysis 
and is stated as follows: 


5(t’) - A sin^rtT 1 + 0^J + 


a 


. /2jt , 

sin fc T + 


(A14) 


where A and 02 are given "by equations (A10) and (All) and the remaining 
parameters are defined as follows : 

_ i - A am 0. I 


Ifl- t 2 


— - A sin 0. . 
0 r l 


- A sin 0-^ - A cos 0 


01 = 2:rrT 2 “ 02 


(A15) 


If it is assumed for the present that and are known, the output 
waveform of the system for the situation involving limiting can he established 
as a piecewise continuous function. For 0 < t % t-^, the initial condition 

must he changed from i to - i and the argument in equation (AlU) must he 
replaced hy 2 jtt ! = 2rtT + - 2jrr 2 ) as is apparent from figure ik. Therefore, 


5(t) = -A sin^2nr + rt - 0 2 ) - B^xp - - Tg + ~j ^ sin | 


sin f ( T - t 2 + i) + 0 . 


(° = T = T i) (ai6) 
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During the next time interval, the system is on the positive limit so that 

S(t) = i ( T 1 ^ T ^ T 2 ) (A17) 

The response during the next interval of time can he obtained by substituting 
2itT ' = 2itT - 2rtT2 into equation ( Al4 ) : 


6(t) 


= A sin ^2 7tr 


*2) 


+ B<exp 





( T 2 = T = t 3) (A 18 ) 

In a similar manner, 

5(t) = - i (Tj % T <; TjJ (A19) 


Finally, for the last interval the sign on the initial condition must again he 
changed to - i and the argument in equation (Al4) must he replaced hy 
2 ^:t t = 2jtt - 



f 

r_ 



r~ — 

l 

t- 

K 

OJ 

p) - B<exp 

r -¥f- 
n v 

- T b) 

> sin 

|( T - T M + 0 o 


l 

L_ 

— 

J 

LX ^ 


where from figure l4 


T, = Tn + ~ 

3 1 2 


T 4 = T 2 + 2 


< t ^ 1.0) (A20) 


(A21) 

(A22) 


The times Tg and may now be determined. From equation (A 6) it is 

apparent that 


notation) applied to the system whereas 


— J sin 2 jtt represents a disturbing force (in nondimens ional 

a/ / 0 \2 

=£■] 5 may be interpreted as a non- 

dimensional elastic restoration or spring force. During limiting, the spring 

force equals j jo. If a > 1.0, that is, |3 0 > the disturbing force is 

greater than the spring force and the system will remain on the limit until the 
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disturbing force equals the spring force. Equating the spring force (fig. 1*0 
to the disturbing forces produces an equation for time Tq. 


t q = sin-l(±) (0 > 1.0) (A23) 

Since the input is a periodic function, the time at which the system departs 
from the positive limit is given by the following equation: 


T 2 “ 2 ~ T 0 


(A24) 


If Tg is known, it is possible to return to equation (Al8) and solve for t ^ 
since 



r 

— — - 

1 

r— — 

B(t) = A sin^jrr^ - 0, 

3 ) + B/exp 

- ^( r 3 " T 2) 

\ sin 

§( T 3 “ t 2) + 00 


(A25) 


If, however, cr < 1.0 (B 0 < (3-,), the disturbing force is always less than 

V “ l > /o_\ 2 t 


or equal to the maximum spring force 


(. 


« 7 * 


It should be noted that for 


this situation it is still possible for the system to limit because of the 
overshoot associated with low damping and/or a near -resonance frequency ratio. 

A cusped response is obtained in this case as illustrated in figure 4(c). Under 
this condition, periodicity requires that 


T 3 



(A26) 


and by using equations (A25) and (A26), both t 2 a 11 '! can be determined. 

The analytical solution was checked by simulating the original constant- 
coefficient differential equations (eqs. (A 6) to (A8)) on an analog computer. 
The results are presented in figure h for several combinations of the input 
parameters. As is apparent from the figures, the comparison is very good. The 
sinusoidal response equations for both the linear and nonlinear modes are sum- 
marized in the section "Derivation of the Describing Function.” 
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APPENDIX B 


NUMERICAL DATA FOR THE MAXIMUM DYNAMIC PRESSURE CONDITION 


The following numerical values were used for the coefficients required in 
equations (l) and (2): 


f S' (x)dx = -8.955(10 ) ~3 
d n 


Tb - D - f s' (x)dx = 1.056(10)3 

J n 


m m 


£ /V 

m J 0 


x) S' (x)dx = -6.92 


T 

— = if-2.7 
m 


Hi f (x - x) S' (x)dx = -9.764(10)“ 6 

I J n 


2 r ^ 

/ (x - x) S^xjdx = -0.1992 

i J n 


£V 

I 


r 1 

/ (x - x) 2 s' (x)dx = -2.655(l0) -2 
J 0 


T„x . 

-j- = 7-29(10) 


-2 
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fr Tbe aeronautical and space activities of the United States shall be 
conducted so as to cotitribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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